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Abstract 

We compute time delays for gravitational lensing in a flat ACDM Swiss cheese universe. We 
assume a primary and secondary pair of light rays are deflected by a single point mass condensation 
described by a Kottler metric (Schwarzschild with A) embedded in an otherwise homogeneous 
cosmology. We find that the cosmological constant's effect on the difference in arrival times is 
non-linear and at most around 0.002% for a large cluster lens; however, we find differences from 
time delays predicted by conventional linear lensing theory that can reach ^ 4% for these large 
lenses. The differences in predicted delay times are due to the failure of conventional lensing to 
incorporate the lensing mass into the mean mass density of the universe. 

PACS numbers: 98.62. Sb, 95.30.Sf, 98.80.-k 
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While investigating a recently proposed effect of the cosmological constant A on gravita- 
tional lensing [1], the authors discovered that angular deflections caused by the embedding 



of point masses in 



as a few percent {2]. The difference in the deflection angles is primarily caused by making 
the deflector's mass part of the universe's mean mass density rather than as an addition 
to it, as is assumed in conventional linear lensing theory. In this letter we use the same 
model to compute analytic expressions for time delays between images produced by these 
lenses. To incorporate the lenses into the mean homogeneous mass density and to enable 
us to evaluate time delays beyond the linear term we use a simple Swiss cheese cosmology 



lomogeneous cosmologies can differ from the Einstein value by as much 



1 . Estimations of these time delays using a similar lensing model have recent^ 



y appeared 



in [6|, 0] and earlier estimates, some using less precise models, appeared in j8l-|l2|. These 
attempts primarily focus on determining the cosmological constant's effect on lensing, in- 
cluding time delays, whereas this analytic work emphasizes the significantly larger effect 
caused by embedding the lens into the cosmology rather than the conventional approach of 
simply superimposing the lensing mass on top of the homogeneous mean density. 

Swiss cheese models are exact general relativistic (GR) solutions which replace comoving 
homogeneous spheres in Friedman-Lemaitre- Robertson- Walker (FLRW) universes by spher- 
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ical condensations (lenses) Because a Swiss cheese lens is part of the mean density 

there is a range beyond which it ceases to deflect passing light rays. That range is just the 
comoving radial boundary Xb of the homogeneous sphere that is replaced by the condensa- 
tion. Beyond that distance the gravity caused by a condensation and a homogeneous sphere 
are the same. 

In this calculation we assume the background FLRW cosmology is a flat {Q^ + Q\ = 1) 
ACDM model. The simplest Swiss cheese cosmology is constructed by replacing the co- 
moving spheres of homogeneous dust by point masses js, 4]. When A 7^ the metric in an 
evacuated void is the Kottler metric IJ, Q (Schwarzschild with a cosmological constant 
present). In \^ we calculated in detail the deflection angle of a photon passing through 
a Kottler condensation. We now compute analytic expressions for the time delay caused 
by encountering such a deflector. We compute the diflFerence in arrival times ATq of two 
light rays emitted at the same time from two sources at equal comoving distances from the 
observer. One ray is assumed to travel entirely in FLRW and arrive at time Tq and the 
other ray is assumed to encounter a deflector and arrive at the observer at the later time 
To = To + ATo (see Figm. 

The Kottler metric 14'1 can be written as 



ds' = --f{r)-'c'dt' + 7(r)'rfr' + r\de' + sin^ 9 1 



where 7 ^(r) = wl — /3^(r), 



/?V) = ^ + ^, 

r 6 



(2) 



and Tg = 2Gm/c^ is the Schwarzschild radius of the condensed mass. The flat FLRW metric 
for the background cosmology can be written as 



-c^dT^ + R{Tf dx^ + x\de^ + sin^ ed4>^) 



with its time development determined by 

R 



R 



= H{R) = Ho\nA + n, 



Ro 

'r 



(3) 



(4) 



To satisfy the junction conditions of GR the boundary of the comoving sphere of dust 
removed [x6 = constant; T, 0, = arbitrary] is matched to an expanding sphere of radial 
timelike geodesies [(^^(T), r5(T)) with 6^,0 arbitrary] in the Kottler vacuole, by requiring 



f^m^(i?oX6)', n{T) = R{T)Xb. 



(5) 



FIG. 1. Swiss cheese gravitational lensing. Bottom: the spatial comoving paths of three photons. 
The deflected photon leaves a source 5, enters a Kottler hole at point 1, exits at 2 at a deflected 
angle a, and then proceeds to the observer at 0. A second photon travels straight from the source 
to the observer as if the Kottler condensation were absent. An imagined third photon also travels 
as if the condensation were absent but is reflected at point B by angle a before arriving at 0. 
The point B (the reflection point) is the intersection of the forward and backward extensions of 
respective FLRW rays SI and 20 drawn as if the Kottler hole were absent, xsi and X20 are 
comoving distances respectively from the source S to entrance point 1 and from exit point 2 to the 
observer, xso and xdo are comoving distances from the source and deflector respectively to the 
observer. The comoving distances from source to reflection point B and from B to the observer 
are respectively XSB = Xsi + XiB and xbo = XB2 + X20- The angular positions (as seen by the 
observer) of the image and source relative to the optical axis are respectively 9j and 9s- Top left: 
a blow up of the two triangles 1B2 and 1D2. The angles 0i,^i,q:, and Acj) are described below 
Eq. (fTTI) and are the same as those computed in Top right: a space time diagram showing the 
difference in arrival times ATg and ATp of three photons originating from source S at time Ts and 
arriving at observer 0. The upper photon encounters a Kottler condensation at time Ti, exits at 
T2, and then arrives at the observer at time Tq = Tq + ATg + ATp. The lower photon arrives at Tq 
after traveling on a straight line entirely in FLRW. The middle photon arrives at time Tq + ATg 
after traveling on two straight segments SB and BO both entirely in FLRW but whose directions 
differ by the angle a. 
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FIG. 2. Fractional differences between Swiss cheese (ATsc) and conventional Schwarzschild (ATsch) 
time delays as functions of source position 6s in units of Einstein ring radius Oe for a lO^^M© (left) 
and a lO^^M© (right) lens. The solid and dashed curves are the fractional time delay differences in 
the potential and the geometrical parts respectively. The dotted and dot-dashed curves are those 
parts of the fractional differences due to the square root term (contains A) in Eqs. (fMll and (fTOll . 

A consequence of Eq. ([5]) for flat FLRW is that the net dust removed equals the condensed 
Schwarzschild mass. 

A photon (see Fig.d]) traveling in a flat ACDM universe starting from a source at time 
Ts and comoving distance xs'o from the observer, moves along a null geodesic of Eq. (|3j), and 
arrives at time Tq where 

r^o cdT _ fRo cdR 

We want to compare this straight path travel time to the time Tq — Ts it takes a photon to 
travel along a deflected path from the source through a Kottler condensation and then to 
the observer. If the entrance time into and exit time from the Kottler hole are respectively 
Ti and T2, the respective comoving distances in FLRW traveled by the deflected photon 
before and after encountering the Kottler hole are 

r^i cdR ^ f^o cdR 

= l.s M' = In. mi- 

Combining these with Eq. (|6]) we have an equation to solve for the difference in travel times 

[To cdT , fR2 cdR 

L ^ = = + ^ + - 
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Ax is just the difference in the comoving distance the deflected photon would have traveled if 
it went straight in FLRW for the actual travel time Tq — Ts and the comoving distance from 
the source to the observer (a travel time of To — Ts if the deflector were not encountered). 
The difference in comoving distances Ax can be factored into the sum of a geometrical part 
and a potential part respectively Ax^ and Axp by adding and subtracting the sum of two 
comoving distances Xib + Xb2 (see Fig.d]), 



The path SBO is that of an imaginary photon that starts from the source S and travels to 
point 1, just as the deflected photon does, but then continues on in FLRW until it is reflected 
at point B by an angle a (same as the Kottler deflection angle) and flnally arrives at the 
observer at Tq + AT^. The deflected photon, if traveling straight, would go an additional 
comoving distance Axp and consequently arrives an additional time AT^ later. Axp is the 
difference in the comoving distance a photon would have traveled in FLRW during the 
Kottler crossing time T2 — Ti minus the comoving distance the imaginary reflected photon 
did travel in crossing the un-evacuated sphere. We identify AT^ and AT^ as the geometrical 
and potential parts of the time delay for purposes of comparison with standard linear lensing 
delays (l6|. A decomposition is not unique but this one is useful and results in Eqs. ( |TOl) and 



To proceed further we must relate comoving distances Xsb + Xbo and Xib + Xb2 to known 
deflector-observer and source-observer distances Xdo and xs'o as well as evaluate the integral 
in Eq. (|9j). We can avoid computing xsb and Xbo to all but the lowest order Xsb ^ Xsq — Xdo 
and Xbo ^ Xdo by using the law of cosines on the triangle SBO. To the accuracy we give 
ATq it suflices to replace the left hand side of Eq. (JSj) by cATq/Rq. The resulting geometrical 
part of the time delay is 




(9) 




(10) 
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where the deflection angle a has been taken from Eq. (32) of [2j. The source-deflector- 
observer distances in Eq. (fTOl) are the usual angular diameter distances. For small impact 
angles 0i this result reduces to the standard lensing result given in Eq. (12) of [16]. For 
larger impact angles the difference is signiflcant, even vanishing when 0i ^ 7r/2 as the 
photon grazes the Kottler void. 

To evaluate the potential part of the time delay in Eq. ([9j) we need to determine Xib+Xb2- 
If we apply the law of sines and cosines to the triangles 1B2 and 1D2 respectively, the result 
is 



The deflected photon's trajectory angles in this expression are shown in Fig.[T] and are the 
same as those used in [2]. They are deflned as follows: In the plane of the deflected photon 
(spherical polar angle = 7r/2) the incoming photon's slope is tan(^i, the photon impacts 
the Kottler void at point 1 with azimuthal angle tt — 0i, the photon exits the void at point 
2 with slope tan(^i + a) and azimuthal angle 0i + A(/). By deflnition the angles a, ^i, and 
A(/) are negative. Equation (fTTI) can be evaluated using values for these angles given in 
as series that depend on the deflector's mass, the cosmological constant, the photon's 
minimum Kottler radius Tq, and the photon's entrance angle (f)i. We evaluate the integral 
in Eq. ([9]) by expanding it as a series in Ai? = R2 — Ri to order {ARy using Eq. (|4j) where 
AR is the change in the radius of the universe that occurs while the photon traverses the 
Kottler condensation. It is related to the change in the Kottler coordinate of the boundary 
Ar = ri){T2) — ^^(Ti) that occurs during the photon's transit by Eq. ([5j), i.e.. 



The change Ar is related to the change in the azimuthal angle Acj) = 02 — 0i by the orbit 
equation (11) of [2] and A(f) is given by Eq. (fT3l) below. In Eq. (13) of j2| we gave an expression 
for A(j) but not to the accuracy needed to evaluate non-linear corrections to the time delay. 
Using the method described in [2] we computed Acj) to the next higher order. 



XiB + Xb2 = 2x6Cos((/)i + A(t)/2) cos(A(/)/2 - ^ - c^/2)/ cos(a/2). 



(11) 



AR = ^ A 



(12) 
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2 log I cot — > sin 01 (^4 — 3 sin^ tan^ 0i 



1 



36-9 sin 0i - 70 sin^ 0i + 9 sin^ - 59 sin^ 



-81 sin^ 



— sec 01 
9 ^ 



tanc 



+ -At' 



(l + Msin^ ^1 - 12sin^0i) 



2 , f 01 I 2 7 
- log <^ cot y ^ sec 01 



(13) 



In this expression see Eq. ([2]), is the expansion velocity of the Kottler boundary as seen 
by a static observer at the instant the photon enters the void. When the above is inserted 
into the integral in Eq. (jHI) and combined with Eq. ffTTl) the potential part of the time delay 
can be approximated by using a from Eq. (32) and from Eq. (18) of [2]. The result is 



2il±^r, 



log I cot ^ I ~ COS 01 ^1 + - COS^ 01 



■ cot 0iCOS 011 



3 ro V XdJ_ 



(14) 



Arriving at Eq. (fT4l) also required converting the redshift Zi of the photon when just entering 
the Kottler void to the slightly smaller redshift Zd of a source located at the deflector's 
distance. To the accuracy needed here 



(1 + zi) = (1 + Zd)[l + /3i cos 01 + + ^iXb/Xd)]- 



(15) 



The result for ATp in Eq. ( IT41) can be compared with the original linear result given in 
Eq. (17) of flG^. The Swiss cheese result reproduces the original result when 0i ^ (but 
7^ 0). The comparison requires two images seen at small 0i's. The potential part of the time 
delay also vanishes as it should when 0i 7r/2, i.e., when the photon misses the Kottler 
hole. The most signiflcant part of the correction is to the linear term and just as with the 
geometrical term the cosmological constant only appears as a part of the expansion velocity 
of the Kottler boundaries [the square root terms in Eqs. (fTOl) and (fT4l) come from /3i]. 

In Fig. [2] we illustrate the use of Eqs. (fTOl) and (fT^ by comparing conventional and Swiss 
cheese time delays for two point mass lens systems with Zd = 0.5, Zg = 1.0 and masses 
m = 10^^ Mq and 10^^ Mq. The cosmological parameters used are = 0.3, = 0.7 
and Hq = 70kms~^ Mpc~^. We graph the fractional diflFerence in time delays between Swiss 
cheese and standard Schwarzschild lensing as functions of source position 9s. For both lenses 



the abscissa 9s varies from to 2.5 9 e where is the classical Einstein ring radius. For the 
10^^ Mq case the total relative correction to the potential part of the time delay is seen to 
be larger than 2% ranging up to 4% at 2.5 9e- The effect of A is to increase the potential 
time delay by as much as 0.002%. The geometrical part of the time delay differs from the 
conventional Schwarzschild result by as much as 0.1% of which a reduction by a factor of 
3 X 10"'' was caused by the A terms in Eq. (fTOl) . For the 10^^ Mq deflector we observe 
a change in the potential part of the time delay of up to 0.4%. The A term inceases the 
potential part of the time delay by as much as 2 x 10"''. The geometrical part of the time 
delay is different by as much as 4 x 10~^ whereas A's effect is a reduction of no more than 
a factor of 3 x 10"^°. 

For currently observed strong lensing by clusters the net time delay corrections are only 
^0.2%. For example a pair of images for a lens like SDSS J1004+4112 at redshifts Zd = 0.68, 
Zg = 1.73, a lens mass of M = 2.3 x 10^^ and a source position 9s = 0.15^^, produces a 
Schwarzschild time delay of ^7.26 yrs with a Swiss cheese increase of ^ 4.9 days. As 
indicated by Fig. [2] future lensing configurations could yield corrections 10 times larger. 
However, exactly how important Swiss cheese lensing is for lens modeling awaits additional 
theoretical work on the embedding of more realistic distributed mass lenses. 

Work on this project was partially supported by NSF grant AST-0707704 and US DOE 
Grant DE-FG02-07ER41517 and B. Chen wishes to thank the University of Oklahoma Foun- 
dation for a fellowship. 



[1] W. Rindler, & M. Ishak, Phys. Rev. D, 76, 043006 (2007). 

[2] R. Kantowski, B. Chen, & X. Dai, Astrophys. J., in press. [arXiv:0909.3308[ 

[3] A. Einstein & E. G. Straus, Rev. Mod. Phys., 17, 120 (1945). 

[4] E. Schiicking, Z. Phys., 137, 595 (1954). 

[5] R. Kantowski, Astrophys. J., 155, 89 (1969). 

[6] K.-E. Boudjemaa, M. Guenouche, and S.R. Zouzou. ^iv:1006.0080"l 

[7] T. Schiicker, arXi v: 1006.32341 

[8] M. Ishak, Phys. Rev. D, 78, 103006 (2008). 

[9] M. Serene, Phys. Rev. D, 77, 043004 (2008). 

9 



[10] T. Schiicker, & N. Zaimen, A&A, 484, 103 (2008). 

[11] M. Sereno, Phys. Rev. Lett., 102, 021301 (2009). 

[12] T. Schiicker, Gen. Relativ. Gravit., 41, 67 (2009) 

[13] R. Kantowski, Astrophys. J., 155, 1023 (1969). 

[14] F. Kottler, Ann. Phys. (Leipzig), 361, 401 (1918). 

[15] C. C. Dyer, & R. C. Roeder, Astrophys. J., 189, 167 (1974). 

[16] J. H. Cooke & R. Kantowski, Astrophys. J., 195, Lll (1975). 



10 



